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WEDGE  THEORY  /  COMPOUND  MATRICES: 
PROPERTIES  AND  APPLICATIONS 


Debra  L.  Boutin,  M.S.* 
Ronald  F.  Gleeson,  Ph.D. 
Robert  M.  Williams,  Ph.D. 


2  August  1996 


ABSTRACT 

The  Navy  utilizes  matrices  to  analyze  radar  signals  to  determine  the  direction  and 
velocity  of  aircraft.  Matrix  analysis  is  also  useful  in  the  sonar  classification  of  submarines. 
One  powerful  tool  for  obtaining  information  about  matrices  is  wedge  theory.  (The  tradi¬ 
tional  terminology  is  “compound  matrix  theory”,  whereas  modern  texts  speak  of  “map¬ 
pings  on  the  exterior  algebra”.)  Wedge  theory  is  a  fundamental  tool  in  multilinear  algebra 
with  important  applications  to  group  representations  and  tensor  analysis.  Current  research 
indicates  that  it  may  also  be  useful  in  analyzing  noisy  data  matrices,  but  this  potential  has 
not  yet  been  fully  explored.  The  purpose  of  this  report  is  to  collect  details  about  wedge 
theory,  in  one  accessible  place,  to  facilitate  future  exploration  of  this  topic.  First,  basic 
properties  of  the  wedge  operation  axe  given  along  with  definitions  and  examples.  Then,  an 
application  to  calculating  the  rank  of  a  matrix  with  noise  is  considered.  Finally,  since  the 
basic  constructions  can  now  be  easily  implemented  on  desktop  computer  algebra  systems, 
the  procedures  for  several  such  packages  axe  illustrated. 


*  Supported  by  National  Defense  Science  and  Engineering  Graduate  Fellowship. 
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I.  INTRODUCTION 


The  Navy  utilizes  matrices  to  analyze  radar  signals  to  determine  the  direction  and 
velocity  of  aircraft.  Matrix  analysis  is  also  useful  in  the  sonar  classification  of  submarines. 
The  role  of  matrices  in  Navy  data  analysis  and  a  tool  for  extracting  information  from  these 
matrices  have  been  discussed  in  an  earlier  report  (NAWCADWAR-96-21-TR)  by  Gleeson, 
Stiller  and  Williams. 

Another  tool  for  squeezing  information  out  of  matrices  is  wedge  theory.  In  the  older 
literature  wedge  theory  is  referred  to  as  compound  matrix  theory.  In  more  modern  texts 
one  speaks  of  mappings  on  the  exterior  algebra.  Wedge  theory  appears  to  have  the  potential 
to  be  quite  useful  in  analyzing  noisy  data  matrices,  but  this  potential  has  not  yet  been 
fully  explored.  The  purpose  of  this  report  is  to  collect  details  about  wedge  theory,  in  one 
accessible  place,  to  facilitate  future  exploration  of  this  topic. 

The  wedge  product  of  a  matrix  is  defined  in  Section  II.  Also,  certain  basic  facts  and 
properties  of  the  wedge  product  are  discussed.  In  addition,  the  eigenvalues  and  character¬ 
istic  equation  for  the  wedge  product,  along  with  a  method  for  efficiently  computing  the 
coefficients  of  the  characteristic  equation,  will  be  considered.  Finally,  we  examine  theorems 
relating  to  different  types  of  compound  determinants. 

In  Section  III  we  create  a  rank  three  4x4  matrix  with  noise  deliberately  added. 
We  then  show  that  the  wedge  products  of  this  matrix  have  a  proportionally  larger  nullity 
than  the  original  matrix.  This  increased  nullity  may  afford  a  method  for  gaining  greater 
sensitivity  in  determining  the  effective  rank  of  the  matrix. 

Tools  such  as  the  wedge  product  have  become  more  feasible  with  the  development  of 
desktop  computer  algebra  software  packages.  We  will  show  in  Appendix  A  how  to  produce 
the  wedge  products  using  Maple,  Mathematica  and  Fermat. 


II.  THEORETICAL  BACKGROUND 

A.  DEFINITIONS 


Given  a  matrix  A,  our  goal  is  to  define  A P(A),  a  matrix  whose  entries  are  the  pxp 
minors  of  A.  This  matrix  has  been  called  the  pth  “wedge”  of  A,  “compound”  of  A  or 
“exterior  product”  of  A,  depending  on  the  literature.  In  this  paper  we  shall  refer  to  AP(A) 
as  the  pth  wedge  of  A  or  as  the  order  p  wedge  of  A. 

First  a  little  necessary  notation. 
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Definition:  Let  A  be  an  n  x  n  matrix  whose  entry  in  the  row  and  j  column  is 
denoted  by  aij. 


/aU 

<*1,2  •  ’  * 

<*l,n  \ 

<*2,1 

<*2,2  *  ’ - 

<*2,n  j 

^  an,  1 

<*n,2  •  *  • 

<*n,n  ^ 

Lexicographically  ordered  p-sets:  Fix  a  number  p  between  1  and  n,  inclusive.  Consider 
all  possible  sets  of  p  distinct  numbers  between  1  and  n.  Order  the  elements  of  each  set 
by  increasing  magnitude.  Order  the  sets  by  increasing  magnitude  of  the  first  elements, 
and  when  the  first  elements  are  equal,  by  increasing  magnitude  of  the  second  elements, 
and  when  the  first  and  second  elements  are  equal,  by  increasing  magnitude  of  the  third 
elements,  etc.  That  is,  order  the  sets  using  a  lexicographic  or  dictionary  order.  Denote 
these  sets  by  Si,  S2, . .  •  ,  S ^  so  that  Si  <  S2  <  •  •  •  < 

Example:  Let  n  =  4  and  p  —  2.  Then  Si  =  {1,2},  S2  =  {1,3},  S3  =  {1,4},  S4  =  {2,3}, 
S5  =  {2,4},  S6  =  {3,4}. 


Definition:  Let  Ajj  be  the  p  x  p  matrix  formed  by  the  intersection,  in  A ,  of  the  rows 
whose  numbers  are  in  the  set  S,'  and  the  columns  whose  numbers  are  in  the  set  Sj .  That 
is,  Aij  is  formed  by  the  entries  a,k,k  of  A  where  h  is  an  element  of  Sj  and  A:  is  an  element 
of  Sj,  maintaining  the  relative  placement  of  the  entries. 

Note:  The  determinant  of  any  p  x  p  submatrix  of  A  is  called  an  order  p  minor  of  A,  or 
alternately,  a  p  x  p  minor  of  A.  This  terminology  will  be  used  throughout  the  paper. 


Example:  Let  A  be  an  abstract  4x4  matrix.  Let  p  =  2.  Then  S2  —  {1, 3}  and  S5  —  {2,4}, 
so  A2  5  has  entries  from  the  intersection  of  rows  1  and  3  with  columns  2  and  4  of  A.  Thus 


Definition:  Let  Ap( A)  be  the  matrix  of  size  (")  x  (np)  whose  (ij )th  entry  is  the  determinant 
of  the  matrix  Aij. 


Example:  Let  A=(  \  ;  -Then, 

\  0  0  0  2/ 
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Basic  Facts  about  Wedge  Products 


Let  A  be  an  n  x  n  matrix. 

1)  The  matrix  AP(A)  has  dimension  (”)  x  (”).  Thus  the  dimensions  of  the  various 
wedge  products  of  A  correspond  to  the  nth  row  of  Pascal’s  Triangle. 

2)  For  theoretical  purposes,  A°(A)  is  defined  to  be  the  lxl  matrix  with  entry  1. 

3)  A1  (A)  =  A  (A)  =  A. 

This  is  clear  since  by  definition,  A1  (A)  has  entries  which  are  determinants  of  the  lxl 
submatrices  of  A.  That  is,  the  entries  of  A(A)  are  the  entries  of  A. 

4)  A"(A)  =  (det(A)). 

This  is  an  immediate  consequence  of  the  definition  of  A"(A).  The  entries  of  A n(A) 
are  the  determinants  of  the  n  x  n  submatrices  of  A.  But  A  is  the  only  n  x  n  submatrix  of 
itself. 
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B.  PROPERTIES  OF  THE  WEDGE  OPERATOR 


The  beauty  and  usefulness  of  the  wedge  operator  is  demonstrated  in  the  following 
properties.  A  good  reference  for  most  of  these  properties  is  Determinants  and  Matrices 
by  A.  C.  Aitken,  [1].  In  particular,  properties  1),  4),  6),  10)  and  11)  can  be  found  there. 
Property  2)  can  be  found  in  Algebra,  Volume  2  by  P.  M.  Cohn,  [2]. 

Theorem;  Let  A  and  B  be  n  x  n  matrices.  Let  A  be  a  real  number.  Let  I  and  I  be 
identity  matrices  of  appropriate  dimensions.  Then: 

Property  1)  A P(I)  =  I'. 

This  can  be  seen  easily  from  the  definition  of  AP(A). 


Property  2)  f\p(AB )  =  A P(A)  Ap  ( B ). 


Example;  Let  A  = 


so  we  can  compute  A3(A2?)  = 


A \B)  = 
equal  to  A3(AB). 


and  B 


2-3  5-6 

0-3  3  6 

0  0-6  12 
-0  0  0-6 


.  Then  AB  = 


Further,  A 3  (A)  = 


Now  we  can  compute  A3  (A)  A3  ( B )  = 


which  is 


Property  3)  AP(AA)  =  Ap  •  AP(A). 


In  the  case  A  =  I  an  easy  computation  yields  AP(A I)  =  A PI'.  If  A  ^  I,  we  can  write 
\A  =  XI  A.  Using  2),  the  wedge  product  becomes  AP(AA)  =  AP(A  I  A)  =  AP(A  I)  Ap  (A)  = 
XpI'  Ap  (A)  =  Ap  Ap  (A). 


Example:  Let  A  = 


A3  (2  A) 


-P 


-16  -8 
-24 
0 

0  0 


(l  1  1 

[0-1  1 
l  0  0  2 
\0  0  0 
-8  —16  \ 
24  48  \ 

48  48  J 

0  -48/ 


I .  Then  we  can  compute  2 A  = 


which  can  easily  be  seen  to  be  23 A3  (A)  =  8 


2 
2 
2 
6 

-2  -1  -1  -2  \ 
0  -3  3  6  \ 

0  0  6  6  J 
0  0  0  -6  / 


Property  4)  (  Ap  (-4))  1  =  AP(A  l). 
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This  is  a  consequence  of  1)  and  2).  Since  AA  1  =  I,  we  have  that  I1  =  A P(I) 
AP(AA“’1)  =  AP(A)  Ap  (A"”1).  Thus  AP(A~X)  is  the  inverse  of  AP(A). 


Example:  With  A  as  above,  A  1  =  I  * 

\o  o  o 

(-2  -1  -] 

o  o  < 

o  o  ( 

I .  Thus  A3(A”1)  is  actually  the  inverse  of  A3(A). 


Z1 

1 

-1 

3  \ 

1 0 

-1 

1 

2 

M 

0 

0 

1 

2 

-* 

\o 

0 

0 

4/ 

/ 

-2 

-1 

-1  "I 

0 

-3 

3  ( 

- 1 

0 

0 

6  < 

\ 

0 

0 

0  -< 

and  A3(A  x)  — 


o  o  o  - 


i  i  —i  1 
'  2  6  6  6 

0  -I 

0  0  I  I 

o  o  o-i 


10  0  0 
0  10  0 
0  0  10 
0  0  0  1 


Property  5)  (  Ap  (A))*  =  A*(A*). 

This  is  an  easy  consequence  of  the  definition  of  A P(A). 


Example:  With  A  as  above,  A1  =  I  x  x  2  0  J  and  A3(At)  = 


10  0  0 
1-10  0 


\i  1  1  3> 

can  easily  be  seen  to  be  the  transpose  of  A3 (A)  = 


-2  -1  -1  -2 
0-336 
0  0  6  6 
0  0  0  -6 


-2  0  0  0 
-1-3  0  0 

-13  6  0 

-2  6  6  -6 


,  which 


Property  6)  If  A  is  a  symmetric  matrix,  then  AP(A)  is  symmetric. 

This  is  an  immediate  consequence  of  5).  If  A  is  symmetric  A  =  A1.  Then  by  5), 
Ap(A)  =  Ap(At)  =  ( Ap  (A))1,  which  means  that  AP(A)  is  symmetric. 


Example:  Let  A  = 
also  symmetric. 


10  12 
0  10  1 
10  10 
2  10  1 


,  a  symmetric  matrix.  Then  A3  (A)  = 


0-202 
-2-4  2  0 

0  2-4-2 

2  0-20 


Property  7)  If  A  is  upper  (lower)  triangular,  then  AP(A)  is  upper  (lower)  triangular. 

Example:  Note  that  in  the  example  for  3)  the  original  matrix  is  upper  triangular  and  so 
is  the  wedge  product.  The  example  for  5)  is  a  lower  triangular  matrix,  A1,  whose  wedge  is 
also  lower  triangular. 


Property  8)  If  A  is  a  diagonal  matrix,  then  AP(A)  is  diagonal. 
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This  is  an  immediate  consequence  of  7).  If  A  is  diagonal  matrix,  then  it  is  both  upper 
and  lower  triangular;  by  7)  so  is  AP(A). 

(1  0  0  0\  / 6  o  o  o \ 

0  0  3  0  )  •  Then  A3  (A)  =  (  °  o  i°2  o  )1S  a^so  diagonal. 

0004/  V00024/ 


Property  9)  The  complex  conjugate  of  AP(A)  is  the  pth  wedge  of  the  complex  conjugate 
ofA. 


Example:  Let  A 


/»+*' 

1 

2+i 

2  i  > 

V 

a=(  : 

1— i 

-1 

1— 2t 

I  is  the 

V  0 

0 

— t 

— i 

/ 

V  0 

0 

0 

i  ' 

/ 

/ 

-2i 

-2i 

4  i 

5+3  i  \ 

A3(A)  = 

0 

0 

2* 

0 

1  -• 
1+* 

T4i 

V 

0 

0 

0 

1-i  / 

Then  A3  (A)  = 


is  the  complex  conjugate  of  A.  We  can  then  compute 


which  is  the  complex  conjugate  of  A3 (A). 


and 


Property  10)  The  Hermitian  conjugate  of  AP(A)  is  the  pth  wedge  of  the  Hermitian  con¬ 
jugate  of  A. 

Since  the  Hermitian  conjugate  of  a  matrix  is  the  complex  conjugate  transposed,  this 
is  a  straightforward  consequence  of  9)  and  5). 


Property  11)  (Sylvester)  det(Ap(A))  =  (det(A))^p 


Example:  In  the  example  for  3),  det(A)  =  —6  while  det(A3(A))  =  —216  =  (— 6)(2). 


Property  12)  As  a  linear  operator  on  a  vector  space  of  dimension  (”),  AP(A)  is  injective 
(surjective)  if  and  only  if  A  is  injective  (surjective)  as  a  linear  operator  on  a  vector  space 
of  dimension  n. 

For  those  who  are  less  familiar  with  injectivity  and  surjectivity,  we  provide  the  fol¬ 
lowing: 

A  linear  operator  A  is  said  to  be  injective  if  for  any  two  vectors  x  and  y,  Ax  =  Ay  only 
if  x  —  y.  That  is,  no  two  distinct  vectors  have  the  same  image  under  the  linear  operator. 
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A  linear  operator  is  said  to  be  surjective  if  for  any  vector  v  there  is  some  vector  u  so 
that  Au  =  v.  A  consequence  of  A  being  surjective  is  that  the  image  of  the  vector  space 
under  the  linear  operator  A  is  again  the  whole  vector  space. 

For  finite  dimensional  vector  spaces  the  notions  of  injectivity,  surjectivity  and  non-zero 
determinant  are  equivalent.  Thus  Property  12)  follows  from  Property  11). 


C.  EIGENVALUES  OF  THE  WEDGE  PRODUCT 

Theorem:  If  Aj, . . .  ,  An  are  the  eigenvalues  of  A,  then  the  eigenvalues  of  AP(A)  are  the 
(”)  distinct  products  of  the  A i  taken  p  at  a  time. 

Proof: 

Given  a  matrix  A  there  is  related  matrix,  called  the  Jordan  form,  which  has  the  same 
eigenvalues  as  A,  with  the  added  advantage  of  being  lower  triangular.  The  Jordan  form 
will  be  extremely  helpful  in  finding  the  eigenvalues  of  AP(A). 

Let  B  be  the  Jordan  Form  of  the  matrix  A.  Then  B  is  similar  to  A ,  which  means 
that  B  =  PAP -1  for  some  matrix  P.  Since  similar  matrices  have  the  same  eigenvalues 
and  since  the  Jordan  form  is  lower  triangular,  B  has  the  eigenvalues  of  A  as  its  diagonal 
entries. 

We  can  compute  A P{B)  =  A P(PAP~1)  =  A P(P)  Ap  (A)  Ap  (P)-1  by  the  properties  of 
the  wedge  operator. 

This  means  that  A P(B)  and  AP(A)  are  also  similar  matrices  and  therefore  have  identi¬ 
cal  eigenvalues.  Thus  we  have  that  the  eigenvalues  of  AP(A)  are  the  eigenvalues  of  A P(R). 

Note  that  since  B  is  a  lower  triangular  matrix,  A P(B)  is  lower  triangular  and  its 
diagonal  entries  are  |Bi,i|,  \B2, l2|,  •  •  •  ,  in  the  notation  of  Section  II.A.  Thus  the 

eigenvalues  of  AP(B)  (and  therefore  of  A P(A))  are  precisely  these  |J9^|. 

Let’s  see  what  one  of  the  |Sj,i|  looks  like. 

Recall  the  p-sets  Si,...  , defined  earlier.  Let  Si  =  {*1,12,...  , ip}.  Then  is 

the  determinant  of  the  p  x  p  matrix  formed  by  entries  of  B  from  rows  . .  , ip  and 

columns  •  ■  ■  ,  ip.  That  is: 


^«1  ,«2 

bh,is 

bi\,ip 

^«2,*2 

^»2,«3 

bh,iP 

Bi,i  = 

^*3, *2 

bia,ia 

bis  ,iP 

kip  ,»2 

bip,is 

bip  ,ip 
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Since  B  is  a  lower  triangular  matrix,  whenever  ij  <  ik  we  know  that  &<,-,»*  =  0.  In 
Bi  i  this  occurs  precisely  when  j  <  k.  Thus  B{tj  is  a  lower  triangular  matrix.  Further, 
since  the  diagonal  entries  of  B  are  the  eigenvalues  of  A  we  know  that  bj.  ^  ,  the  ij 

eigenvalue  of  A.  Thus, 


Now  we  can  easily  compute  =  A^  Aj2  •  •  •  A ip ,  an  eigenvalue  of  A P(B)  and  therefore 
of  AP(A). 

Since  the  eigenvalues  of  AP(A)  are  precisely  the  which  are  the  products  of  the 

eigenvalues  of  A  taken  p  at  a  time,  the  eigenvalues  of  AP(A)  are  the  (”)  distinct  products 
of  the  A,-  taken  pat  a  time. 

Example:  If  A  is  a  4  x  4  matrix  with  eigenvalues  Ai,  A2,  A3,  A4  then  A2 (A)  has  eigenvalues 
AiA2,  A1A3,  AxA4,  A2A3,  A2A4,  and  A3A4,  and  A3(A)  has  eigenvalues 
AiA2A3,  AiA2A4,  AiA3A4  and  A2A3A4. 

Example:  Using  the  numeric  example  from  Section  II. A,  we  see  that  A  has  —3, —1,1 
and  2  as  its  eigenvalues,  which  implies  that  A 2(A)  has  eigenvalues  —3  •  —1  =  3,  —3  •  1  = 
—3,  — 3  •  2  =  —6,  — 1  •  1  =  —1,  —  1  •  2  =  —2,  and  1-2  =  2.  We  can  see  immediately  that 
this  is  true  by  looking  at  the  upper  triangular  matrix  A2  (A)  which  we  computed. 


D.  THE  CHARACTERISTIC  POLYNOMIAL  OF  THE  WEDGE  PRODUCT 

We  would  like  to  be  able  to  find  the  characteristic  coefficients  (that  is,  the  coefficients 
of  the  characteristic  polynomial)  of  AP(A),  directly  from  the  characteristic  coefficients  of 
A.  This  turns  out  to  be  a  straightforward  task  for  two  of  the  characteristic  coefficients* 

For  the  purposes  of  this  paper,  the  characteristic  polynomial  of  the  matrix  A  is  defined 
to  be  det(a;I  -  A),  though  in  some  literature  it  is  defined  as  det(A  -  xl).  The  polynomials 
resulting  from  these  two  definitions  differ  by  a  power  of  -1.  However,  we  find  it  convenient 
to  have  the  highest  degree  term  in  the  characteristic  polynomial  be  positive,  and  our 
definition  ensures  that. 

Let  c*  be  the  coefficient  of  the  term  of  degree  n  -  k  in  the  characteristic  polynomial 
of  A.  In  the  characteristic  polynomial  of  Ap( A),  let  dk  be  the  coefficient  of  the  term  of 
degree  (J)  -  k. 
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More  explicitly,  the  characteristic  polynomial  of  A  will  be  written  as  xn  +  c\  xn  1  + 
•  •  •  +  cn-\x  +  cn,  while  if  we  let  m  =  (”),  the  characteristic  polynomial  of  AP(A)  will  be 

written  a s  xm  +  dixm~l  -) —  •  +  dm-\x  +  dm. 


It  is  a  well  known  fact  in  linear  algebra  that  in  the  characteristic  polynomial  of  an 
n  x  n  matrix,  the  coefficient  of  the  term  of  degree  n  —  k  is  (— 1  )kSk  where  s*  is  the  kth 
symmetric  function  on  the  eigenvalues  of  the  matrix.  That  is,  sjt  is  the  sum  of  products 
of  eigenvalues  taken  k  at  a  time. 


So  in  the  characteristic  polynomial  of  AP(A),  d\,  the  coefficient  of  the  term  of  degree 
(”)  —  1,  is  (— l)1^!,  the  negative  of  the  sum  of  the  eigenvalues  of  AP(A).  Since  these 
eigenvalues  are  products  of  eigenvalues  of  A  taken  p  at  a  time,  d\  is,  up  to  sign,  actually 
the  pth  symmetric  function  on  the  eigenvalues  of  A.  Thus: 


Theorem:  In  the  characteristic  polynomial  of  AP(A),  d\  =  (— l)p+1cp. 


We  now  use  a  similar  method  to  find  the  constant  term,  dyiy 

Recall  from  the  properties  of  the  wedge  product  that  det(Ap(A))  =  (det(A))^p-1^  and 
that  the  constant  term  of  the  characteristic  polynomial  of  a  matrix  is,  up  to  sign,  the 
determinant  of  the  matrix.  Thus: 


Theorem:  In  the  characteristic  polynomial  of  AP(A),  the  constant  term,  d^y  is 
up  to  sign. 


The  other  characteristic  coefficients  of  AP(A)  don’t  lend  themselves  to  such  simple 
formulae.  However,  they  are  still  polynomials  in  the  cjt’s,  and  we  can  actually  find  them. 


Method  for  Finding  Characteristic  Coefficients 


Let  A  be  an  abstract  n  x  n  matrix  with  characteristic  polynomial 


p(x)  =  xn  +  cizn  1  H - hCn-!X  +  C  „. 


Suppose  that  B  is  a  matrix  with  the  same  characteristic  polynomial  as  A.  Then  the 
eigenvalues  of  B ,  i.e.  the  roots  of  the  characteristic  polynomial,  are  identical  to  the  eigen¬ 
values  of  A.  Thus  Ap(A)  and  AP(B)  also  have  identical  eigenvalues.  Since  the  coefficients 
of  the  characteristic  polynomial  of  a  matrix  are  completely  determined  by  its  eigenvalues, 
A P(B)  has  the  same  characteristic  polynomial  as  AP(A). 

So  what  we  want  is  a  matrix  B ,  with  the  same  characteristic  polynomial  as  A ,  whose 
entries  are  0’s,  l’s  and  cjfc’s.  If  AP(B)  also  has  entries  of  0’s,  l’s  and  ct’s,  the  characteristic 
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coefficients  of  A P(B)  will  be  polynomials  in  the  c*’s.  Luckily,  the  companion  matrix  for 
p(x)  fits  all  these  requirements  perfectly. 


(° 

0  ••• 

0 

1 

0  ••• 

0 

Cn— 1 

Let  B  = 

0 

1  ••• 

0 

—  Cn- 2 

Ko 

• 

0  ••• 

1 

-Cl  / 

the  companion  matrix  for  p(x). 


B  has  the  same  characteristic  polynomial  as  A,  and  the  entries  of  B  and  A P(B)  are 
0’s,  l’s  and  c*’s,  so  the  characteristic  coefficients  of  A P(B)  (and  therefore  of  AP(A))  are 
polynomials  in  the  c*’ s. 

In  particular  we  can  use  a  computer  algebra  program  such  as  Mathematica,  Maple  or 
Fermat  to  compute  the  characteristic  polynomial  of  A P(B).  The  result  will  yield  a  formula 
for  each  characteristic  coefficient  of  A P(A)  in  terms  of  the  coefficients  of  the  characteristic 
polynomial  of  the  original  matrix,  the  c*’s. 

This  needs  only  be  done  once  for  each  n  and  each  p.  Then,  given  a  matrix  with  a 
known  characteristic  polynomial  we  need  only  use  the  appropriate  formulae  to  find  the 
characteristic  polynomial  of  the  p ^  wedge.  Some  of  these  formulae  are  listed  in  Appendix 
B. 

Example:  Let  A  be  a  4  x  4  abstract  matrix  which  has  characteristic  polynomial  of  the 

r  /  0  0  0  —  C4  \ 

form  x4  +  cix3  +  c2x2  +  c3x  +  c4.  If  we  let  B  =  I  J  J  J  I J  I ,  we  can  then  compute 

\0  0  1  -Cl  ) 


/° 

0 

c4 

0 

0 

°v 

0 

0 

0 

0 

0 

c4 

0 

/° 

— c4 

0 

1 

0 

0 

0 

-c2 

0 

0 

0 

c3 

c4 

0 

and  A 3(B)  = 

0 

0 

0 

0 

1 

0  X? 

0 

c4  1 

0 

1 

-Cl 

0 

0 

C3  I 

Vi 

-ci 

C2 

c3  J 

Vo 

0 

0 

1 

-ci 

C2/ 

We  can  now  compute  the  characteristic  polynomial  of  A2 (I?),  which  is  x6  +  (  c2)x  + 
(c3Cl  - c4)x4  +  (2c4c2  -  c4cj  -  4)x3  +  (c4c3C!  - cl)x2  +  (-c\c2)x  +  cl ,  and  the  characteristic 
polynomial  of  A3 (13),  which  is  x4  +  c3x3  +  c4c2x2  +  c2Cix  +  c3. 

Thus  we  have  the  characteristic  polynomials  of  A2(A)  and  A3 (A)  respectively,  and  the 
characteristic  coefficients  are  polynomials  in  the  characteristic  coefficients  of  A. 

Note  that  in  the  characteristic  polynomial  of  A2(A),  d\  =  (— l)2+1c2  and  d§  =  (c4)^ 
as  promised. 
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Also,  in  the  characteristic  polynomial  of  A3(A),  c?i  =  (— 1)3+1C3  and  d4  =  (c4)(*), 
again  as  expected. 

Example:  Let  A  be  the  4x4  numeric  matrix  from  the  example  in  Section  II. A. 

The  characteristic  polynomial  of  A  is  x4  +  x3  -  7x 2  -  x  +  6,  so  cj  =  1,  c2  =  -7, 
C3  =  —1,  and  c4  =  6. 

Using  these  in  the  formulae  above,  we  get  that  the  characteristic  polynomial  of  A2 (A) 
is 


x6  +  -(-7)x5  +  ((-1X1)  -  (6))x<  +  (2(6)(— 7)  -  (6)(1)2  -  (-l)V 
+  «6)(— 1)(1)  -  (6)2)x2  +  (-(6)2(-7))x  +  (6)3 
=  x6+  7x5  -  7x4  -  91x3  -  42x2  +  252x  +  216, 

and  the  characteristic  polynomial  of  A3  (A)  is 

x4  +  (-l)x3  +  (6)(— 7)x2  +  (6)2(l)x  +  (6)3  =  x4  -  x3  -  42x2  +  36x  +  216. 


E.  DIFFERENT  TYPES  OF  COMPOUND  DETERMINANTS 


In  this  section,  we  cover  some  results  about  minors  of  “compound  matrices.”  The 
wedge  product  is  one  such  compound  matrix,  but  there  are  other  compound  matrix  con¬ 
structions  different  from  the  wedge  product.  An  excellent  reference  for  this  section  is, 
Determinants  and  Matrices  by  A.  C.  Aitken,  [1]. 

Adjugate  Wedge  Products 

Recall  that  the  entries  of  AP(A)  are  the  order  p  minors  of  A.  Here  we  wish  to  study 
a  matrix,  which  we  will  call  the  pth  adjugate  wedge  of  A,  whose  entries  are,  up  to  sign, 
order  n  —  p  minors  of  A.  This  new  matrix  is  not  quite  A n~p(A),  but  it’s  close.  Many  of  the 
signs  of  the  entries  will  be  different,  and  their  positions  within  the  matrix  will  be  different. 
Specifically,  we  wish  the  entries  of  this  new  matrix  to  be  the  cofactors  of  the  minors  which 
comprise  AP(A).  So  first  we  must  define  what  we  mean  by  the  cofactor  of  a  minor. 

Definition:  Let  A  be  an  abstract  n  x  n  matrix.  Let  M  be  the  minor  which  is  the 
determinant  of  the  mxm  submatrix  obtained  from  rows  r\ , . . .  ,  rm  and  columns  C\ , . . .  ,  cm . 
Then  the  complementary  minor  to  M  is  the  determinant  of  the  (n  —  m)  x  (n  —  m)  matrix 
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obtained  by  deleting  rows  r\ , . . .  ,rm  and  columns  c\ , . . .  , cm.  The  cofactor  of  the  minor 
M  is  the  complementary  minor  times  (— l)ri+'"+rm+Cl+  +Cm. 


Example:  Let  A 


a\  d2  <*>3 
b\  62  63  64 
Cl  C2  C3  C4 
d 1  ^2  ds  d4 


M  is  j1  j4  and  the  cofactor  of  M  is  (— 1) 

«1  «4 


Let  M  = 


_ 1  \2+3+l+3 


|.  Then  the  minor  complementary  to 

61  64  _ /  -1  \  \  t>i  64 

A,  A,  —  V  dA  * 


Definition  Let  the  pth  adjugate  wedge  of  a  matrix  A,  denoted  adjp(A),  be  the  (”)  X  (”) 
matrix  whose  ( i,j)th  entry  is  the  cofactor  of  the  order  p  minor  which  is  the  (j,i)th  entry 
of  AP(A). 


Example:  Let  A  — 


Then 


a3(A)  = 


fli  a2  fl3  <*4 
b  1  b 2  b  3  64 
Cl  C2  C3  C4 

d 1  d%  ^3  d4 


ai  a 2  as 
61  62  ^3 

Cl  C2  C3 

ai  a2  03 
b 1  62  &3 
rfl  <^2  ^3 
a\  a 2  03  j 

Cl  C2  C3 

di  ^2  dz 

b 1  62  63 
Cl  C2  C3 
di  <£2  dz 


adj3(A)  = 


ai 

a2 

04 

ai 

a3  04 

02 

03 

04 

61 

62 

64 

61 

63  64 

62 

&3 

64 

Cl 

C2 

c4 

ci 

C3  c4 

c2 

C3 

C4 

ai 

a2 

a4 

ai 

03  04 

02 

03 

O4 

bi 

62 

64 

61 

63  64 

^2 

63 

64 

di 

d7 

d4 

di 

dz  d4 

d7 

<*3 

d4 

ai 

a2 

a4 

ai 

03  04 

02 

O3 

a4 

Cl 

C2 

C4  1 

Cl 

Cz  C4 

C2 

C3 

C4 

di 

d2 

d4 

di 

dz  d4 

d7 

^3  ^4 

hi 

62 

64 1 

61 

63  64 

62 

&3 

64 

Cl 

C2 

C4 

Cl 

c3  c4 

C2 

C3 

C4 

d\  d<2  d4  I 

di 

^3 

d7 

<*3 

d4 

/ 

<*4  - 

C4 

64  —04  \ 

.  f 

—dz 

C3  ‘ 

-63  03  1 

* 

l 

di  — 

C2 

b7  —0*2  J 

\ 

-dx 

Ci  - 

-61  ai  / 

,  and 


Note  that  adjp(A)  is  “almost”  An_p(A).  Up  to  a  factor  of  (-1),+-7,  the  entry 

of  adjp(A)  is  the  (n  -  j  +  l,n  -  i  +  l)th  entry  of  A n~p(A).  That  is,  to  get  adjp(A),  we 
take  An-p(A),  reflect  the  entries  over  the  right-leaning  diagonal,  and  adjust  the  entries  by 


Example:  Let  A  = 


4  2-2  1-5-2 

0  -2  -8  -1  -4  2 

0  0  2  0  1  0 

0  0  0  -2  -8  -5 

0  0  0  0  2  1 

0  0  0  0  0  -1 


2  10  2 
0  2  1-1 
0  0-1-4 

0  0  0  1 


and  A3  (A) 


Then  we  can  compute  A1  (A)  =  A,  A2(A)  — 


-4  -16  -10  -9 
0  4  2  1 

0  0-2-1 

00  0-2 
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Then  adj*(A)  = 


adj2(A)  = 


(-2119 

o  -2  -2  -10 

0  0  4  16 

0  0  0  -4 

-1  -1  -5  0  2  2  \ 

0  2  8  1  4  -5 

0  o-2  0-1-1  • 

0  0  0  2  8  -2 
0  0  0  0  -2  -2 
000004/ 

1  4  -1  -2\ 

0-1-101.  ••11 


0  0  0  2 


has  the  appropriate  relationship  with  A3 (A); 


is  related  in  the  desired  way  to  A4  2  (A)  =  A2 (A);  and 


adj3(A)  =  I  o  Q  2_i  I  Is  similarly  related  to  A1  (A)  =  A 


Note:  The  first  adjugate  wedge  of  A,  adj*(A)  is  just  called  the  adjugate  of  A  and  is  often 
denoted  simply  by  adj(A).  In  some  literature  this  is  called  the  classical  adjoint  of  A. 

Theorem:  (Cauchy)  det(adj(A))  =  (det(A))"  1. 

Example:  In  the  example  above,  we  can  easily  see  that  det(A)  =  —4  while  det(adj(A))  = 
—64  =  (— 4)4-1  which  fits  our  theorem. 

Theorem:  A-1  =  a™ ^adj(A). 

Example:  With  A  as  in  the  example  above 


dsW^>  =  -; 


-2119 
0  -2  -2  -10 
0  0  4  16 

0  0  0  -4 


0  0 
0  0 


-1  -4 


det(A) 


adj(A)  A 


-1  -4 


0  2 
1  -1 
-1  -4 
0  1 


10  0  0 
0  10  0 
0  0  10 
0  0  0  1 


Theorem:  (Jacobi)  Let  A  be  a  matrix.  Then  any  minor  of  order  r  of  adj(A)  is  equal  to 
the  cofactor  of  the  corresponding  minor  in  the  transpose  of  A,  multiplied  by  det(A)r_1 . 

Example:  In  the  previous  numeric  example,  we  can  see  that  the  minor  of  adj(A)  of 
order  3,  obtained  by  deleting  the  last  row  and  last  column,  is  16.  The  cofactor  of  the 
corresponding  minor  of  the  transpose  of  A  is  (_ij1+2+3+1+2+3  .1  =  1.  When  this  is 
multiplied  by  det(A)3-1  =  (— 4)2  we  get  back  the  value  of  the  original  minor. 
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Theorem:  (Franke)  Any  minor  of  order  r  in  the  pth  wedge  of  A  is  equal  to  the  co¬ 
factor  of  corresponding  minor  in  the  transpose  of  the  pih  adjugate  wedge  of  A,  times 

(det(A))r  ^  p  \ 

Example:  Again  let  A  = 

(10  0  0 

: : 

-2  0-1  2 

Consider  the  order  2  minor  of  A3  (A)  which  is  the  determinant  of  the  matrix  obtained 
from  rows  1  and  2  and  columns  2  and  3.  It  is  ^  **  =  8.  The  cofactor  of  the 

corresponding  minor  in  (adj3(A))t  is  (— 1)1+2+2+3  °  =  (— 1)2+3+1+2  •  (— 2)  =  —2.  We 

can  see  that  det(A)  is  —4  and  8  =  (—2)  •  (— ■ 4)2  (»),  as  predicted. 

The  Bazin  Hybrid 

There  are  more  ways  of  creating  new  matrices  from  old  ones  which  result  in  relation¬ 
ships  similar  to  those  we  saw  above.  In  the  following,  we  discuss  Bazin  compound  matrices 
and  their  more  general  form,  Reiss  compound  matrices,  and  theorems  regarding  both. 

Definition:  Let  A  and  B  be  n  x  n  matrices.  The  Bazin  hybrid  compound  of  A  and 
B  is  the  n  X  n  matrix  whose  ( i,j)th  entry  is  the  determinant  of  the  matrix  obtained  by 
replacing  column  i  of  A  with  column  j  of  B. 

Example:  Let  A  and  B  be  abstract  2x2  matrices  with  general  entries  aij  and  bij, 
respectively.  Then  the  Bazin  Hybrid  of  A  and  B  is 

b\,i  ai  ,2  &1,2  ^1,2 

&2,1  0*2,2  &2,2  <22,2 

0tl,l  &1,1  <*1,1  &1,2 

<22,1  &2,1  012,1  &2,2 

Note:  There  is  no  reason  to  believe  that  Bazin  hybrid  of  A  and  B  is  equal  to  the  Bazin 
hybrid  of  B  and  A.  In  fact  they  are  not  generally  equal.  So  we  call  the  Bazin  hybrid  of  B 
and  A  the  dual  Bazin  hybrid  of  A  and  B. 

Theorem:  (Bazin)  The  determinant  of  the  Bazin  hybrid  compound  of  A  and  B  is  equal 
to  (det(A))"  1det(B). 
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Example:  Let  A 


of  A  and  B 


to  (det(A))n  1det(B)  =  (6) 


,4-1 


0  \  /2  0  0  0  \ 

J  ),£=  (  °  ;  •  Then  the  Bazin  hybrid  compound 

2/  \0  1  1  -1/ 

,  which  has  determinant  864.  We  can  see  that  this  is  equal 
•  (4),  as  was  stated  above. 


The  dual  Bazin  hybrid  (that  is  the  Bazin  hybrid  of  B  and  A)  is 


with 


determinant  384  =  (det(B))”  1det(A)  =  44  1 
reverse  the  roles  of  A  and  B. 


6.  This  also  fits  Bazin’s  theorem  when  we 


Theorem:  (Reiss)  Any  minor  of  order  r  of  the  Bazin  hybrid  compound  of  A  and  B 
is  equal  to  the  cofactor  of  the  corresponding  minor  in  the  transpose  of  the  dual  hybrid, 
multiplied  by  (det( A)) r  1(det(2?))r  "+1. 


Example:  Let  A  and  B  be  as  above.  Then  again  the  Bazin  hybrid  compound  of  A  and 

12  0  0  ox  /2  —2  0-6< 

®  6  4  o  )  ^he  transpose  of  the  dual  Bazin  hybrid  is  I  ®  *  Jj  2 

-12  -3  1-3/  Vo  0  0-8- 


B  is 


is 


The  minor  of  the  Bazin  hybrid  obtained  from  rows  3  and  4  and  columns  1  and  3 
=  54.  The  cofactor  of  the  corresponding  minor  of  the  transpose  of  the  dual 


6  4 
-12  1 


-2 

-4 


hybrid  is  obtained  from  rows  1  and  2  and  columns  2  and  4  and  is  (— 1)3+4+1+3 

—1  •  (—36)  =  36.  We  can  compute  (det(A))r  1(det(B))r  "+1  to  be  (6)x(4)-1  =  |  and 
multiplying  36  by  this  gives  us  54,  the  value  of  the  original  minor. 


The  Reiss  Hybrid 


Note:  This  is  a  simple  generalization  of  the  Bazin  hybrid.  Similarly,  the  theorems  axe 
generalizations  of  those  for  the  Bazin  hybrid. 

Definition:  Let  A  and  B  be  n  x  n  matrices.  Define  the  pth  Reiss  hybrid  compound 
of  A  and  B  to  be  the  (”)  x  (")  matrix  whose  (i,j)th  entry  is  the  determinant  of  the 
matrix  obtained  by  replacing  columns  *i, . . .  ,ip  of  A  with  columns  j\,. . .  ,jp  of  B,  where 
Si  =  ip}  and  Sj  =  (ji , . . .  ,jp}  are  the  ordered  p- subsets  of  {1,2,...  ,n}  defined 

at  the  beginning  of  this  paper. 

Note:  We  can  also  define  the  pth  dual  Reiss  hybrid  of  A  and  B  as  the  pth  Reiss  hybrid 
compound  of  B  and  A. 


Example:  Let  n  =  4  and  p  =  2.  Let  A  and  B  be  abstract  matrices  with  general  entries  a{j 
and  bitj.  Then  Si  =  {1,2},  52  =  {1,3},  S3  =  {1,4},  S4  =  {2,3},  S5  =  {2,4},  56  =  {3,4}. 
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The  entry  in  row  1 ,  column  5  of  the  second  Reiss  hybrid  of  A  and  B  is  the  determinant  of 
the  matrix  with  columns  1  and  2  of  A  replaced  with  columns  2  and  4  of  B : 


&1,2 

&1,4 

<*1,3 

Ol,4 

&2,2 

t>4,4 

^2,3 

02,4 

&2,3 

h,4 

03,3 

03,4 

&4,2 

&4,4 

04,3 

04,4 

Theorem:  (Reiss,  Picquet)  The  determinant  of  the  pth  Reiss  hybrid  compound  of  A 
and  B  is  equal  to  (det(A))^  r  ^(det(B))^p_l^. 

Notice  that  when  p  =  1  we  have  the  statement  of  the  theorem  by  Bazin. 


r 


Example:  Let  A  and  B  be  as  in  the  previous  example.  Then  the  second  Reiss 

-12  0  0  0  0  O' 

12  8  o  o  o  o 

-6  2  -6  0  0  0 

0  -4  0  -4  0  0 

-9  -1 

9  9  -3  3  -3 

(3\ 


hybrid  compound  is 


:\ 


3-1  3  0  . 

-3  -2/ 


As  stated  in  the  theorem,  it’s  determinant 


is  13824  =  (det(A))®  (det(£))®  =  63  •  43. 


The  dual  Reiss  hybrid  compound  is 


/ 


\ 


-2  0  0  0  0 

3  3  0  0  0 

3  1  -4  0  0  0 

-3  -3  0  -6  0  0 

-9-14-2  8  0 

9  9  0  -6  -12  -12 


°\ 

0  ' 


It’s  determinant  is 


also  13824  =  (det (B))®  (det(A))^  =  43  :  63. 

Theorem:  (Reiss,  Picquet)  Any  r  x  r  minor  of  the  pth  Reiss  hybrid  compound  is 
equal  to  the  cofactor  of  the  corresponding  minor  in  the  transpose  of  the  dual  Reiss  hybrid, 

multiplied  by  (det(A))r  ^p-1^(det(B))r  ^  r  \ 

Example:  Let  A  and  B  be  as  before.  Then  we  need  the  second  Reiss  hybrid  of  A  and  B 

-12  0  0  0  0  0\  /  -2  3  3  -3  -9  9  \ 

/  a  a-ioi  n  1 


/— 12  0  u  u  u  u\ 

'  12  80000' 
-6  2  -6  0  0  0 
0  -4  0  -4  0  0 

-9-1  3-1  3  0 

9  9  -3  3  -3  -2/ 


,  and  the  transpose  of  the  dual  Reiss  hybrid 


\ 


0  3  1-3-1  9 

0  0  -4  0  4  0 

0  0  0  -6  -2  -6 

0  0  0  0  8  -12 

0  0  0  0  0  -12 


Consider  the  minor  of  the  second  Reiss  hybrid  which  is  the  determinant  of  the  matrix 
obtained  from  rows  and  columns  1,2,3  and  4.  It’s  value  is  —2304.  The  cofactor  of  the 
corresponding  minor  in  the  transpose  of  the  dual  is  (— 1)1+2+3+4+1+2+3+4(— 96)  =  —96. 

Computing  (det(A))4-^~^(det(R))4-('  2  ^  yields  (6)1  •  (4)1  =  24  and  multiplying  -96  by 
this  gives  us  —2304,  as  expected. 
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III.  EXAMPLE  OF  USEFULNESS 


In  the  earlier  report  (NAWCADWAR-96-21-TR)  by  Gleeson,  Stiller  and  Williams  it 
was  shown  that  the  characteristic  coefficients  (c*’s)  could,  be  used  to  predict  the  effective 
rank  of  a  noisy  matrix.  These  coefficients,  when  properly  normalized,  fall  below  predeter¬ 
mined  threshold  values  for  k  greater  than  the  effective  rank.  These  normalized  coefficients 
are  called  P*’s. 


To  illustrate  how  the  wedge  product  could  be  used  in  this  type  of  analysis,  we  create 
the  following  matrix: 

28.05  6.45  10.33  \ 

-11.56  -15.03  -36.72  | 

33.83  37.56  38.39  I  ' 

-37.40  -17.85  -19.93/ 


17.91 
-5.97 
22.04 
V -24.17 


The  matrix  A  was  generated  by  first  creating  a  rank  three  4x4  matrix,  and  then 
adding  a  small  random  noise  contribution  to  each  element.  For  readability,  each  element 
was  then  rounded  to  two  decimal  places.  The  details  of  the  generation  process  are  spelled 
out  at  some  length  in  NAWCADWAR-96-21-TR. 

For  the  matrix  A ,  we  have: 

Pi  =  0.53,  P2  =  0.60,  P3  =  0.57,  P4  =  0.09. 


The  earlier  report  studied  7x7  matrices  with  different  effective  ranks  and  noise  levels. 
For  small  noise,  the  threshold  was  found  to  be  typically  in  the  range  of  0.2  and  0.3.  Let  us 
assume  that  this  threshold  range  does  not  vary  drastically  and  applies  to  4  x  4  matrices. 
The  above  distribution  of  the  Pt’s  has  P3  above  the  threshold  and  P4  below  the  threshold. 
This  is  the  profile  of  a  matrix  whose  effective  rank  is  three. 

Now  let  us  consider  the  second  wedge  of  A.  The  second  wedge  of  A  is  the  following 
6x6  matrix: 


-39.57 

-230.62 

-595.90 

-347.01 

910.63 

-81.47  \ 

-12.43 

530.56 

459.88 

835.58 

727.55 

-140.47 

8.23 

-163.86 

-107.32 

-259.65 

-172.86 

55.88 

52.85 

107.08 

580.25 

74.30 

798.60 

802.25 

-56.15 

-256.63 

-768.42 

-355.64 

-1142.77 

-355.87 

-6.64 

514.42 

488.52 

800.76 

761.42 

-63.36  / 
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Recall,  that  if  the  given  matrix  A  has  rank  three,  then  only  three  of  the  four  eigenvalues 
are  nonzero.  With  noise  added  the  fourth  would  also  be  nonzero,  but  small  compared  to 
the  other  three.  The  eigenvalues  of  the  second  wedge  are  equal  to  the  products  of  A’s 
eigenvalues  taken  two  at  a  time,  so  there  are  (*)  or  three  significant  eigenvalues  in  A2(A). 
With  three  significant  eigenvalues,  A2  (A)  should  appear  to  have  effective  rank  three. 

The  Pjt’s  for  the  A2  (A)  are  the  following: 

Pa  =  0.60,  P2  =  0.56,  P3  =  0.57,  P4  =  0.22,  P5  =  0.13,  P6  =  0.09. 


Assuming  that  the  7  x  7  thresholds  work  for  6  x  6  matrices,  P5  and  P6  clearly  fall 
below  the  threshold  range.  P4  is  on  the  border;  whereas,  P3  is  higher  than  the  threshold 
range.  This  profile  indicates  that  A2 (A)  is  either  rank  three  or  four.  Actually,  the  fact 
that  P6  is  lower  than  the  threshold,  alone  implies  that  the  rank  of  A  is  very  likely  less  than 
four.  Moreover,  when  we  see  that  P5  is  also  less  than  the  threshold,  the  likelihood  of  the 
rank  of  A  being  less  than  four  is  amplified.  Finally,  P3  is  greater  than  the  threshold.  This 
implies  that  the  rank  of  A2 (A)  is  at  least  three.  This,  in  turn,  implies  that  the  rank  of  A 
itself  is  three.  The  main  point  here  is  that  A2 (A)  has  more  coefficients  with  which  to  work 
and  has  greater  nullity.  Using  the  normalized  characteristic  coefficients  of  A2  (A)  should 
give  us  a  greater  handle  and  enhanced  sensitivity  in  determining  the  rank  of  the  matrix  A. 

As  discussed  in  Section  II,  it  is  not  necessary  to  actually  compute  A 2(A)  to  determine 
its  Pit’s.  All  one  needs  are  the  characteristic  coefficients  of  A2 (A)  ,  and  these  can  be 
computed  in  terms  of  polynomials  in  the  characteristic  coefficients  of  A.  A  few  examples 
of  these  polynomials  are  included  in  Appendix  B. 
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IV.  CONCLUSION 


This  report  may  be  regarded  as  a  primer  on  the  theory  of  wedge  products. 

In  Section  II  we  have  brought  together  definitions,  properties  and  theorems  relating 
to  the  wedge  product.  Also,  we  have  discussed  how  the  companion  matrix  can  be  used  to 
compute  the  coefficients  of  the  characteristic  polynomials  of  the  various  wedge  products. 
Finally,  theorems  which  relate  to  other  forms  of  compound  matrices  have  been  explained. 

In  Section  III  we  have  illustrated  how  the  wedge  product  has  the  potential  to  be  useful 
in  the  determination  of  the  effective  rank  of  a  noisy  data  matrix. 

Future  work  in  this  area  includes  determining  quantitatively  the  extent  of  this  po¬ 
tential  for  predicting  the  effective  rank.  That  is,  we  should  determine  whether  or  not 
the  normalized  coefficients  of  the  characteristic  polynomial  of  A2  (A)  are  more  successful 
than  those  of  A  in  finding  effective  rank.  Characteristic  coefficients  for  higher  order  wedge 
products  should  also  be  studied.  Ultimately,  the  order(s)  of  the  wedge  product(s)  that 
optimize  results  and  computation  speed  should  be  compared  with  existing  methods  for 
predicting  effective  rank. 
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APPENDIX  A 

SAMPLE  RUN  SESSIONS 


Here  we  have  included  sample  code  that  illustrates  the  syntax  for  evaluating  wedges 
of  matrices  on  three  different  mathematical  software  packages:  Maple  (version  5,  release 
3.0),  Mathematica  (version  2.2),  and  Fermat  (version  1.1).  In  the  symbolic  examples  below 
WA2  is  the  second  order  wedge  of  matrix  A. 

Maple 

In  Maple,  there  is  a  wedge  product  function  (&')  in  the  differential  forms  library. 
Unfortunately,  the  intent  of  this  function  was  for  differential  forms,  not  for  linear  algebra. 
The  wedge  product  command  does  not  work  on  all  matrices;  specifically,  it  may  fail  on 
matrices  containing  zero  entries.  This  would  put  limitations  on  the  number  of  matrices  for 
which  the  Navy  could  find  the  wedge  product.  In  light  of  this,  George  Nakos  has  written 
a  program  using  Maple  syntax  that  computes  the  wedge  product  of  any  matrix.  This  file 
can  be  written  in  any  text  editor.  Maple  can  access  this  file  with  the  command  ‘read  < 
filename  >;  ’.  The  program  is  as  follows: 

with(linalg) : 
with(combinat ,cartprod) : 

#This  function  takes  a  list  of  lists  [[a,  b, ...],  [al,61, ...], ... 

#and  computes  the  Cartesian  product  [a,  6,...]  x  [al,  61, ...]  x  ... 
lesslistsO  :=  proc(lis) 
local  11, car; 

11  :=  □; 

car  : =  carprod(lis) ; 
while  not  car [finished]  do 

11  :=  [op(ll) ,car [nextvalue] ()]  od: 

RETURN (11); 
end: 

#  This  function  takes  a  list  of  numbers  and  returns  a  list  of  lists 
#of  numbers  with  entries  <  to  the 

# corresponding  entries  of  the  original.  For  example, 

#Lesslists([2, 4, 3]);  yields  [1, 1, 1],  [1, 1, 2], ...,  [2,4, 3] 

Lesslists  :=  proc(lis) 
local  ll,i,n; 

11  :=  []; 

for  i  from  1  to  nops(lis)  do 
11  :=  [op(ll),[$l..lis[i]]]: 
od: 
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11  :=  lesslistsO(ll) ; 

RETURN (11); 
end: 

#This  function  tests  whether  the  length  of  f  is  n 
lengthNQ  :=  proc(f,n)  evalb(nops(f )  =  n)  end: 

#  From  a  list  lis,  pick  the  elements  of  length  n 
PickLengthN  :=  proc(lis,n)  select  (lengthNQ, lis, n)  end: 

#This  function  deletes  repeated  elements  in  a  list 
DeleteRepeated  :=  proc(lis) 
local  11, i; 

for  i  from  2  to  nops(lis)  do 

if  evalb(not  member(lis[i]  ,11))  then 
11  :=  [op(ll) ,lis[i]] ;fi: 
od: 

RETURN (11); 
end: 

#  This  function  takes  two  >0  integers  n  and  k,  n  >  k,  and  returns  all  pairs 

#of  the  form  (i,j)  with  i  <  j  and  i  <  k  and  j  <  n. 

Submatrixlndex  :=  proc(n,k) 
local  tt; 

tt :=  [$(n-k+l) . .n] ; 

LessLists(tt) ; 
map (convert, ",  set) ; 

PickLengthN (" ,k) ; 

DeleteRepeated(") ; 

tt  :=  map (convert,", list) ; 

Return (tt); 
end: 

$  This  function  takes  a  list  of  numbers  and  returns  the 
^complete  Cartesian  product  of  the  element  lists 
List Of Pairs  :=  proc(lis) 
local  ll,nn,i,j; 
nn  :=  nops(lis); 
for  i  from  1  to  nn  do 
for  j  from  1  to  nn  do 

11  :=  [op(ll) , [lis[i] ,lis[j]]] ; 
od: 

od: 
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RETURN (11); 
end: 

#  MAIN  FUNCTION.  Computes  the  k  x  k  minors  of  matrix  A 
Minors  :=  proc(A,k) 

local  n, dimmat, kk, 11, i; 
n  :=  vectdim(rov(A,l)) ; 

Submatrixlndex(n.k) ; 
dimmat  :=  nops("); 
kk  :=  ListOf Pairs ("") ; 

11  □  ; 

for  i  from  1  to  nops(kk)  do 

11  :=  [op(ll) , det ( submatrix (A, op (kk[i] )))] ; 
od: 

RETURN (matrix (dimmat, dimmat, 11)) ; 
end: 

Example: 


To  compute  the  second  order  wedge  for  the  matrix  A  = 


’  al  a2  a3  a4 ' 
61  62  63  64 
cl  c2  c3  c4 
dl  d2  d3  d4  ■ 


read  in  the  above  file  and  then  type: 

>  A  :=  matrix^ 3, 3,  [al,  a2,  a3,  a4, 61, 62, 63, 64,  cl,  c2,  c3,  c4,  dl,  d2,  d3,  d4]); 

>  W  A2  :=  Minor s(A,  2); 

The  output  to  this  would  be  the  following  matrix: 


,  one  would 


WA2  := 


/al62-a26l  al63-a36l  al64-a461  a263-a362  a264-a462  a364-a463  \ 
alc2— a2cl  alc3— a3cl  alc4— a4cl  a2c3— a3c2  a2c4— a4c2  a3c4  — a4c3  ' 
ald2—a2dl  ald3-a3dl  ald4-a4dl  a2d3-a3d2  a2d4-a4d2  a3d4-a4d3 
61c2— 62cl  61c3— 63cl  6lc4-64cl  62c3-63c2  62c4-64c2  63c4-64c3 
6ld2— 62rfl  61d3— 63dl  61rf4-64dl  b2d3-b3d2  b2d4-b4d2  b3d4-b4d3 
\c\d2-c2d\  cl<*3— c3dl  cld4-c4dl  c2d3-c3d2  c2d4-c4d2  c3d4-c4d3  / 


Mathematica 

Mathematica  has  a  built  in  command  for  computing  the  wedge  product,  namely 
“Minors[A,p]’\  This  command  takes  two  parameters;  the  first  is  the  name  of  the  matrix, 
and  the  second  is  a  positive  integer  denoting  the  order  of  the  wedge. 

Example  : 
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As  in  the  previous  example,  if  we  wish  to  compute  the  second  wedge  of  a  4  x  4  matrix 
in  Mathematica,  we  would  enter: 

A  —  {{al,  0.2,  a3,  a4},  {a5,  a6,  a7,  a8},  {a9,  alO,  all,  al2},  {al3,  al4,  al5,  al6}  } 

WA2=  Minors[A,  2] 


Fermat 

Fermat  is  a  mathematical  software  system  written  by  Robert  Lewis.  While  there  is 
no  built-in  function  for  computing  the  wedge  product,  there  will  be  soon.  Michael  Hirsch 
has  written  a  program  in  the  Fermat  language  that  can  be  used  in  the  absence  of  a  wedge 
product  function.  Run  times  axe  faster  than  Mathematica  and  Maple  for  large  symbolic 
matrices.  Once  the  code  is  within  the  Fermat  shell,  run  times  using  Fermat  should  be  even 
shorter  than  they  are  currently  (using  the  Hirsch  program). 

;This  is  the  main  function  for  the  program 
: Wedge (p,d, mat rix2;  m,n,i)  = 

:cols  =  Cols[p]; 

:rows  =  Deg [p] /cols; 

:n  =  Q  (cols.d); 

:m  =  Q  (rows,d) ; 

:ri  =  0; 

: ci  =  0; 

:b[m,n]; 

:  x  [d]  ; 

:  y  Cd]  ; 

:temp[d,d] ; 

Genr(rows,d,l,l) ; 

:  [b]  =  ™  [b]; 

(D  [temp]  ; 

: [matrix2]  =  [b] ; 

®[b]  ; 

0.; 

;This  function  generates  all  of  the  possible  row  combinations 
;for  the  wedge 
:Genr(m,dd, j ,s;i)= 
if  (dd  =  0) 

then  (:ri+;  :ci  =0;  Genc(cols ,d, 1 , 1)  ) 
else  (  for  (  :i=s,m+l~dd)  do 
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(  :x[j]  =  i;  - 
Genr(m,  dd-1,  j+1,  i+1  ); 

) 

).; 

;This  function  generates  all  of  the  possible 
; column  combinations  for  the  wedge 
:Genc(n,dc,k,l;i)= 
if  (dc  =  0) 

then  (:ci+;  Dump(ri,ci)  ) 
else  (  for  (  :i=l,n+l-dc)  do 

(  :y[k]  ■  i; 

Genc(n,  dc-1,  k+1,  i+1  ); 

) 

)  • ; 

;This  function  calculates  the  determinant  of  each 
; submatrix  for  the  entries  of  the  wedge 
: Dump (row, col;  q,w)s 
:temp[d,d]  ; 
for(  :q=l,  d  )  do 

(  for(  :w=l,d  ))  do 

(  :temp[q,w]  =  p[x[q],y[w]]  ) 

(  ); 

:b  [col,  row]  =  Det[temp]; 

®  [temp]  . ; 

After  reading  this  program  into  Fermat,  using  the  wedge  function  is  easy.  Just  type  the 
following: 

>  Wedge(matrixl  ,d,matrix2) 

where  matrixl  is  the  original  matrix,  d  is  the  order  of  the  wedge,  and  matrix2  is  the 
variable  that  will  be  assigned  to  the  wedge. 

Example: 

In  Fermat  to  compute  the  second  wedge  of  a  4  x  4,  we  would  enter: 

>:  <*[4,4] 

>:  [a]  =  [[al,a2,  a3,a4,  „ 
a5,a6,a7,  a8,  „ 
a9,  alO,  all, al2,  „ 
al3,al4,al5,al6]] 
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>:  u>a2[6,6] 

>  Wedge([a], 2,  [wa2]) 
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APPENDIX  B: 

CHARACTERISTIC  COEFFICIENT  RELATIONSHIPS 

This  Appendix  contains  a  few  sample  coefficients  of  the  characteristic  polynomial  of 
the  wedge  products  (djt’s)  expressed  as  polynomials  in  the  coefficients  of  the  characteristic 
polynomial  (c*’s)  of  the  given  matrix  (A). 

If  A  is  a  3  x  3  matrix: 

For  A2 (A): 

d\  =  —Ci 
d.2  —  C1C3 
dz  =  -c\ 

For  A3 (A): 
d\  =  C3 

If  A  is  a  4  x  4  matrix: 

For  A2 (A): 
d\  —  — C2 

C?2  =  C1C3  —  C4 
dz  =  2C2C4  —  CjC4  —  c2 
C?4  =  C1C3C4  —  c\ 
d5  =  -c2c2 
d6  —  c\ 
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4x4  matrix  cont. 

For  A3  (A): 
d\  —  C3 
=  C2C4 
dz  =  cic^ 
d4  =  c\ 

For  A4  (A): 

d\  =  —c4 

If  A  is  a  5  x  5  matrix: 

For  A2  (A): 
d\  =  — C2 
C?2  =  clc3  ”  C4 

^3  =  C1C5  +  2c2c4  -  c\c4  -  c\ 
d4  =  C3C5  —  3c!C2c5  +  c3  C5  —  C4  +  C1C3C4 
d5  =  — c2c2  -  c2  c3  C5  +  2c2c3c5  +  2C1C4C5  -  2c| 
d6  =  c\  +  ci c2c4c5  -  3c3c4c5  -  +  c2c^ 

=  c4c2  +  2cic3c2  -  c\c\  -  cic|c5 
d6  =  c2c4cj  -  cic3 

dg  =  —  C3C5 

dio  -  c\ 
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5x5  matrix  cont. 

For  A3  (A): 
d\—cz 

d2  =  c2c4  —  C1C5 

<*3  =  —  C4C5  -  2C1C3C5  +  c\c5  +  C\c\ 

d4  =  C2C3  -  cjcj  -  3c3c4c5  +  C1C2C4C5  +  cf 
=  2c\  -  2C1C4C5  -  2c2c3c§  +  c?c3c§  +  C2C4C5 
^6  =  Cic3c4c|  -  cjcj  +  c\c\  -  3ciC2C3  +  c3c| 

di  =  C3C3  +  Cj  C4C5  —  2C2C4C3  —  C1C5 

<*8  =  CiC3c|  -  c4c| 
dg  =  C2c\ 

d\o  =  c| 

For  A4  (A): 
d\  =  — C4 
c?2  =  C3C5 

dz  =  ~c2c\ 
d4  =  ci  cl 
d5  =  c4 
For  A5  (A): 
d\  =  c5 
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If  A  is  a  6  x  6  matrix;- 
For  A2 (A); 
d\  =  —c2 

d2  —  C\C%  —  C4 

d3  =  —  c\  —  C2C4  +  2C2C4  +  C1C5  —  C6 

d4  =  C1C3C4  -  c\  +  c\c5  -  3cic2c5  +  C3C5  -  c\c6  +  2c2c6 
C?5  =  —  C2C2  —  C1C3C5  +  2c2C3  C5  +  2C1C4C5  —  2c2 

— CjC6  +  4c2C2C6  -  2c2C6  -  3C1C3C6  +  2c4c6 

C?6  =  C4  +  CiC2C4C5  —  3C3C4C5  —  c2c2  +  c2c2 

+cfc3C6  —  3ciC2C3C6  +  3c2Cg  —  C2C4C6  +  3C!C5C6  —  2c| 

d7  =  — C\ C4C5  -  +  2cic3c|  +  C4c|  -  cfc2c4c6 

+2c2C4C6  +  C1C3C4C6  —  C4C6  +  cfc5C6  —  CiC2C5C6 

-3c3c5c6  -  c\c\ 

dg  =  c2c4c\  -  cjCg  +  cfc^ce  -  2c2c|c6  +  cic5jc5c6 

— 2c2C3C5C6  —  C2  C3  C5  C6  +  C1C4C5C6  +  C5C6  —  C2C2Cg 

+c^c§  +  3cic3c| 

C?9  =  —  C3C5  —  C1C2C4C5  C<>  +  3C4C5C6  +  C2C2C6  +  C2CjC6 

— c^Cg  +  3cic2C3c|  -  3c|c|  -  c5c4c|  -  3cic5c|  +  2c| 

d10  =  c|  +  C1C3C5C6  -  4c4C2c6  +  c\c4c\  -  2cic3c4cg  +  2 c\c\ 
— 2cjc2csCg  +  3c3C5Cg  +  2c2Cg  —  2c2Cg 

dll  =  — C1C5C6  —  C2C3C5  Cg  +  3ciC4C5Cg  +  CgCg  +  C^Cg 

cic3Cg  -  2c4Cg 

d\2  =  C2c\cl  +  C2Cg  -  2c2c4c|  -  dC5c|  +  c| 
dl3  =  C2c|  -  C3C5C6 
di4  =  C4Cg 
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6x6  matrix  cont. 

For  A2  (A): 

dis  =  —c\ 

For  A3  (A): 

d\  =  C3 

C?2  =  C2C4  —  C1C5  +  C6 

<^3  C4  "4"  ^2^5  2cxC3C5  C4C5  "t"  ^03^6 

^4  =  C3  +  C1C2C4C5  -  3C3C4C5  -  c\c\  +  C2C2  +  C^C6 
3C1C2C3 C6  “I-  3CjCg  "j~  Cj C4Cg  C2C4^6  I  CjCgCg 

<*5  =  c2cjc5  +  C2c3c|  -  2c2c3c2  -  2C1C4C5  +  2c3 

+CxC2C4C6  -  2 cfc3c4c6  -  C2C3C4C6  +  3cic2c6  -  2c2c2c5c6 
+3c^c5c6  +  2cic3c5c6  -  4c4c5c6  +  2c3  c2  -  4cic2c| 

+3c3Cg 

d6  =  cic3c4c§  -  c4cl  +  c\cl  -  3cic2c3  +  czc\ 

+C2C2C6  -  2C1C3C4C6  +  2c3C6  +  C2C2C3C5C6  -  2c2C3C5C6 
— C1C3C5C6  —  3CJC4C5C6  +  6ciC2C4C5C6  —  2C3C4C5C6  +  4c2C2C6 
-c\c\c\  +  2 c\c\  +  c3c3Cg  -  2cic2c3Cg  +  3  c2  c| 

+4cfc4Cg  —  8c2c4c2  -  4cic5c|  +  3cg 

d7  =  c2c3  +  c\c4c\  -  2c2c4c3  -  Cjc|  +  cic2c3c4c5c6 

-3 c|c4c5c6  -  3c2c2c5c6  +  5 c2c|c5c6  +  cfc2c|c6  -  3cic|c|c6 
C2 C3C2C6  +  C2c3c|c6  +  4c  j  C4  C2  Cg  +  c|c6  +  cfc|c| 

— 3cic2c2cg  +  3c3c2  —  2cfc2c4c2  +  5cic2c4c2  +  c2c3c4c2 
-5c2c3c4Cg  +  cjc2Cg  -  cfc5c|  +  4c2c2c5c2 
+C2C5C6  —  cic3c5Cg  —  6c4c5Cg  +  c3Cg  -  6cic2Cg 
+3c3Cg 


31 


NAWCADPAX— 96-220-TR 

6x6  matrix  cont. 

For  A3(A): 

ds  =  cic3c|  -  c4c|  4-  c2c§c§c6  +  c?c2c4cijc6  -  2clc4clc6 

— 4c4 c3 c4 Cg c$  +  4c|c|c6  -  2cic2c3c6  +  cfc|c4c|  -  2c2c\c4c\ 

-2 c\c2c\c\  +  4 444  -  c4Cg  +  c4c3c5c|  -  4 cfc2c3c5c| 
+4cic|c5c^  -  2c3c4c5Cg  +  llcic2c4c5c|  -  5c3c4c5c|  +  2444 
+c2c%cl  -  C j C2 Cg  +  4CjC2Cg  -  44  -  5ciC2C3Cg 
+3C3C3  +  c^c4Cg  -  7c2c4c|  -  5cic5cg  +  3c| 

dQ  =  C2C\  +  CiC2C34c6  ~  6c2C4c|c6  -  Cic|c6  +  C1C3C5C6 

+cfc3c4c5c^  -  5cic2c3c4c5c|  -  c|c4c5c|  -  cfc^c5c|  +  10c2C4C5c| 

— cjc^CgCg  -  2c\cz44  +  5c2c3c§cg  +  6cic444  +  44 
-ClC2cl4  +  44  +  cfc4c|  -  6cfc2C4Cg 
+10cic|c4cg  +  5cjC3c4Cg  -  10c2c3c4Cg  -  bci44 
— CjC5Cg  +  6cfc2C5Cg  —  5C2C5C3  —  3ciC3C5Cg  —  5c4c5c| 

+cf  4  —  5cic2cg  +  6c3c£ 

did  =  4  +  c2c5c6  -  6c4c|c6  +  cfcfci?c;?  -  Ac\c4c\4 

— 2cic3c4cgc§  + 10  444  -  2c!c244  +  ^44  +  44 

-4c2C3C4Cg  +  444  -  4cjC2c4Cg  +  6444  +  4cic3c4Cg 

-4 44  ~  2cfc2c3C5Cg  +  4 c|c3c5c^  -  2cic|c5cj?  -  2 cfc4c5cj? 
+12c4c2c4C5Cg  -  12c3c4c5cg  +  2c4C5Cg  -  4c2c|cg  +  c\4 
— 6c4c20g  +  lOcjCjCg  —  4c2c|  +  6c3c3c|  —  12cic2c3cg 
+6c3Cg  4c4c4c6  2c2c4c6  4"  ScjCsCg 
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6x6  matrix  cont. 

For  A3 (A): 

dn  =  c2cjjc6  +  cic2c3c|c|  -  6c2c4CgCg  -  ci 44  +  ciclc5cl 

+c5c3c4c5c|  -  5cic2c3c4c5c|  -  cfc4c5Cg  -  c^c|c5c|  +  10c2c|c5Cg 
“C1C2C5C6  -  2CjC3C^C^  +  5c2C3CgCg  +  6c!C4CgCg 
+clc3  -  CiC2CgCg  +  C3Cg  +  CjC4Cg  -  6c3C2C4Cg 
+  10ciC2C4Cg  +  5CjC3C4Cg  —  10c2C3C4Cg  —  5CiC4Cg  —  cfc5cj 
+6CjC2C5Cg  -  5C2C5Cg  -  3ciC3C5Cg  —  5c4C5Cg  +  C3Cg 
-5cic2Cg  +  6c3Cg 

d  12  =  Cic3c|c|  -  cAc\c\  +  C2<?zc\c\  +  CjC2 c4 c|cg  -  2c\c4clc\ 

4cic3c4c5c6  H“  4c4c5c6  2cxc2CgCg  4-  c4 c3 c4Cg  2c2CjC4Cg 

— 2CjC2C4c|  +  4c2c|cg  -  C4Cg  +  CjC3C5Cg  -  4CjC2C3C5Cg 
+4cic§c5Cg  -  2 c3c4c5c|  +  llcic2c4c5Cg  -  5c3c4c5cg  +  2cf  cfcg 
+c2clcj  -  c\c2c\  +  4c4C2Cg  -  c\c\  -  5cjC2C3Cg 
+3c|cg  +  c\c4c\  -  7c2c4c®  -  5cic5Cg  +  3cJ? 

rfi3  =  444  +  c\c4c\cl  -  2c2c4c|cg  -  cic\c\  +  cic2c3c4c5c| 

— 3c|c4c5Cg  -  3c4c4c5cg  +  5c2c4c5cg  +  c\c244  -  3cic^ci?c| 
-clc3cl4  +  c2czc\c\  +  4C!C4CgCg  +  44  +  444 
-3cxc2c|c|  +  3c^cj?  -  2c3c2c4cg  +  5cicfc4cf  +  c?c3c4c| 

— 5c2c3c4cg  +  Ci c4 Cg  -  4c$cl  +  4cfc2c5c|  +  c|c5c|  ~  ci czc$4 
6c4c3Cg  -|-  CjCg  6cxc2cg  -f*  3c3Cg 

d\  4  =  Cl  C3C4CgCg  -  444  +  c3c3c|  -  3CiC2CgCg  +  C3  Cg  Cg 
+444  -  2cic3c4Cg  +  2c4c|  +  cfc2c3c5cjj  -  2c|c3c5c| 

-cic|c5c|  -  3c3c4c5cg  +  6c1c2c4c5Cg  -  2c3c4c5cg  +  4c2CgCg 

-444  +  2c2c6  +  cfc3C^  -  2ciC2C3Cg  +  3c|cg 
+4cfc4cg  -  8c2c4cg  -  4C1C5C®  +  3 c\ 
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6x6  matrix  cont. 

For  A3 (A): 

dis  —  C2C^CsCg  +  C1C3C5CI  —  2C2C3C5CQ  —  2C1C4C5C6  +  2c5c6 
+C1C2C4C3  —  2CjC3C4Cg  —  C2C3  C4Cg  +  3c  1  44  —  2c  j  C2  C5  Cg 
+3c|c5cJ?  +  2C1C3C5C6  -  4c4C5Cg  +  2c3Cg  -  4cic2Cg 

+3c3Cg 

C?16  =  c\c\  +  C1C2C4C5CQ  —  3c3C405Cg  —  c\c\c\  +  C2c|cg 

+44  —  3oiC2C3Cg  +  3c|c  g  +  CjC4Cg  —  C2C4Cg 

+ClC5cl 

dn  =  CiC^Cg  +  C^C5Cg  -  2cic3c5c^  -  C4C5Cg  -  C1C2C® 

+3C3C3 

dis  =  c2c4cf  -  cic5c|  +  Cg 
dig  =  C3C® 

d20  =  Cg° 

.  For  A4(A): 

d\  =  — C4 

d2  =  +(0305)  —  C2C6 

ds  =  —c2c\  -  C3C6  +  2c2c4c6  +  C1C5C6  -  Cg 

d4  =  +(cic|)  +  C2C3C5C6  -  3cic4c5c6  -  c§c6  -  44 
+C1C3CI?  +  2c4c| 

d5  =  — c4  -  Ci C3C5C6  +  4c4CgC6  -  4c*cl  +  2cic3c4cg 

-2 44  +  2cic2c5c|  -  3c3c5Cg  -  2cfc|  +  2c2cj! 
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6x6  matrix  cont. 

For  A4  (A): 

C?6  =  +(c3c|c6)  +  c1c2c4c5c%  -  3 c3c4c5c|  -  c\c\c\  -  c2c\cl 
+4 Cg  %CiC2CZc\  -f-  3Cg Cg  -f-  CjC4Cg  +  3ci C5 Cg  2cg 

d7  =  -c2c4c\cl  +  CxCgC^  -  cfcjcg  +  2c2c|4  -  ci c^c5c| 
+2CjC3C5Cg  +  C2C3C5Cg  -  Cj C4C5Cg  -  CgCg  +  cfc2c| 

c2c|  —  3cjc3c| 

d&  =  +(cic^c5Cg)  +  cjcjcl  -  2cic3c|cg  -  c4c\c\  +  c\c2c4c\ 

2c2C4Cg  CiC3C4Cg  “|“  C4Cg  CjCsCg  "f  Ci C2C5Cg 
+3C3C5CI  +  CjCg 

^9  =  -c4c|  -  Cic2c4c5c|  +  3c3c4c5c|  +  c^c§c|  -  c2c§c| 

— Cj c3 Cg  +  3cic2c3c|  -  3c|c|  +  cfc4cjj  -  3cic5c| 

+2c® 

^10  =  +(c2c4Cg)  +  CjC3c5Cg  -  2c2c3c5Cg  -  2cic4c5Cg  +  2c|cg 
+c\c%  4c  j  c2  Cg  +  2c2Cg  +  3cic3c|  —  2c4Cg 

dll  =  C1C3 c4 Cg  +  c4Cg  -  cJc5Cg  +  3cic2c5Cg  -  c3c5Cg 

+44  -  2c2Cg 

dl2  =  +(C3CD  +  C?C4c|  -  2c2C4Cg  -  ClC5Cg  +  c| 

C?i3  CiC3Cg  "f  C4Cg 

du  =  +(c2c|) 

di5  =  -4° 

For  A5  (A): 

d\  =  c5 

d2  =  C4C6 

C?3  —  C3Cg 
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6x6  matrix  cont. 

For  A5 (A): 

d\  —  C2Cg 
d$  —  cic6 
d6  = 

For  A6 (A): 
d\  —  C5 
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APPENDIX  C 
dk  DEPENDENCE 


As  mentioned  in  Section  II  D,  one  can  find  the  dk's  by  looking  at  the  characteristic 
polynomial  of  the  pth  wedge  of  the  companion  matrix  for  a  given  polynomial.  Appendix 
B  specifies  how  the  c*’s  and  dk's  axe  related.  But  how  are  the  dk  s  related  to  each  other? 
Although  the  c*’s  are  algebraically  independent,  the  dk's  axe  not.  Below  is  code  from 
Mathematica  which  uses  the  Grobner  basis  tool  to  calculate  the  dk  relations.  The  Grobner 
basis  tool  was  applied  to  the  equations  from  Appendix  B  using  the  second  wedge  of  a  4  x  4 
matrix.  Recall  from  Appendix  B,  that  the  second  order  wedge  of  the  4x4  matrix  has  six 
d^s.  In  the  ideal  generated  by  this  tool,  there  are  five  polynomials  that  only  involve  the 
six  dk's.  By  setting  these  polynomials  equal  to  zero  and  simplifying  we  get  five  useful  dk 
relations. 


GroebnerB  asis  [  {  d  1  +  c2, 

d2  +  c4  —  cl  *  c3, 

dZ  +  c4  *  cV2  +  c3"2  —  2c2  *  c4, 

dA  +  c4“2  —  cl  *  c3  *  c4, 

d5  +  c2*  c4~2, 

d6  -  c4~3}, 

{cl,  c2,  c3,  c4,  dl,  d2,  dZ,  dA,  d5,  d6}] 

The  output  to  this  command  is  a  Grobner  basis.  The  first  five  polynomials  in  the 
basis,  however,  depend  only  on  the  dk's.  They  are  as  follows: 

d23d6  —  dA3, 
dldA2  -  d22d5, 
dld2d6  —  dAdb, 
dl2dAd6  -  d2d52, 
dl3d62  —  d53. 

Setting  each  of  these  equations  equal  to  zero  and  reducing  them  provides  a  way  of  examin¬ 
ing  the  relations  between  the  dk's.  In  fact,  for  this  example,  it  turns  out  that  the  equations 
reduce  to  the  following: 

d5  =  dl  *  d&3!3 
dA  =  d2*  d62/3. 

Hence,  this  gives  us  an  easier  way  to  calculate  the  dk's. 
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